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ABSTRACT 

In this paper we model the observed absorption measure distribution (AMD) in Mrk 509, which spans 
three orders of magnitude in ionization level with a single-zone absorber in pressure equilibrium. AMD 
is usually constructed from observations of narrow absorption lines in radio-quiet active galaxies with 
warm absorbers. We study the properties of the warm absorber in Mrk 509 using recently published 
broad-band spectral energy distribution observed with different instruments. This spectrum is an 
input in radiative transfer computations with full photoionisation treatment using titan code. We 
show that the simplest way to fully reproduce the shape of AMD is to assume that the warm absorber 
is a single zone under constant total pressure. With this assumption we found theoretical AMD which 
matches the observed one determined on the basis of 600 ks RGS XMM-Newton spectrum of Mrk 509. 
The softness of the source spectrum and the important role of the free-free emission breaks the usual 
degeneracy in the ionization state calculations, and the explicit dependence of the depths of AMD 
dips on density open a new path to density diagnostic for the warm absorber. In Mrk 509 the implied 
density is of the order of 10 8 cm -3 . 

Keywords: instabilities - radiative transfer - method: numerical - galaxies: active - galaxies: individual: 
Mrk 509 


1. INTRODUCTION 

There is a general consensus based on high reso¬ 
lution X-ray data that the majority of the Seyfert 
galaxies contain ionised absorbing gas in their line 
of sight. After the advancements in the observing 
X-ray instrumentation as Chandra, XMM-Newton, 
Suzaku , the extensive studies of the Warm Absorber 
(hereafter WA) in AGN were performed. Many narrow 
absorption lines from highly ionised elements, detected 
with the use of gratings, provided a g reat oppor- 
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Some physical parameters of WA were estimated, but 
some are still difficult to constrain. From energy shifts 
of line centroids, it was found that the absorbing matter 
is systematically outflowing with velocities in the rang e 
10 2 —10 3 kms -1 (jKaspi et al.ll200lHKaastra et al.ll2002l) . 
Directly from observations, ionic column densities of par¬ 
ticular ions could be found and they are typically of the 
order of 10 15-18 cm -2 . In the next step, by assuming sep¬ 
arate ionisation zone for each ion, it is possible to calcu¬ 
late the equivalent hydrogen column densities, knowing 
ion’s ionisation fraction from photoionisation codes such 
as CLOUDY, XSTAR, SLAB, XABS. Such equivalent hydro- 
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;en column densities are in the range of 10 1 
Steenbrugge et all 120051 ICostantini et al.l[2007t . More¬ 
over, those densities correspond to the continuous change 
of ionization parameter £ ( see Eq.[I]for definition), which 
spans the range of a few decades, log(£) ~ — 1 up to 4. 

Since then, several attempts have been made 
to model continuous ionisation structure of the 
WA in sever al AGN (iH olczer. Behar fe Kaspil 120071 : 
Behai! [20091: IDetmers et al. 1201111 . Furthermore, 


Holczer. Behar fe Kaspil (|2007fl proposed to describe ion¬ 


isation structure of the wind by showing absorption mea¬ 
sure distribution (AMD) obtained by a derivative for¬ 
mula (Eq. 5 in this paper). These authors, for the first 
time, have shown that in the case of the source IRAS 
13349+2438, AMD obtained from observations has deep 
minimum in column density that is consistent with the 
negligible absorption of gas with log(£) between 0.8 and 
1.8. Such deep minima are p resent in AMD of other ob¬ 
jects as well (see for ins tance IBeh ar 2009: IDetmers et all 
120111 [Stern et al.ll2014f) . and they are interpreted as the 
observational evidence for thermal instability in a given 
ionisation and the temperature regime. 

To reproduce the observed AMD theoretically, we 
should consider continuous ionisation structure of WA 
in such a way that we can obtain proper normalisation 
of the AMD func tion. It has been done recently by 
iStern et all (I2014H assuming a radiation pressure con¬ 
finement (hereafter RPC) of the WA material. The au¬ 
thors considered that the WA consists of the stratified 
clouds, for which P gas increases gradually with decrease 
of -Prad- The exponential decrease of the radiation pres¬ 
sure is caused via absorption on the material when the 
distance from the illuminating source, i.e. active nucleus, 
gradually increases. The authors were successful in re¬ 
producing the normalisation and the slope of AMD for 
Syl galaxies by showing the comparison of the model 
with observational points of AMD for the six best stud¬ 
ied sources. Nevertheless, they were not able to quanti- 
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tatively reproduce the deep minimum in column density 
for those six objects present in the log£ between 0.8 and 

2 . 

In this paper, we show how AMD for the Syl.5 galaxy 
Mrk 509 can be successfully reproduced by the WA be¬ 
ing under the condition of constant total pressure, w hich 
is a sum of gas and radiatio n pressure (iRozahsk a et al.l 
120061 : iGoncalves et alJf2006fl . The radiation pressure is 
computed from the radiation field at each point of the 
computations and it is added to the gas pressure. We 
neglect other pressure components, for instance the mag¬ 
netic one. We used the broad band spectral energy distri¬ 
bution of Mrk 509 from the paper of lKaastra et al.l (120111 ) 
(hereafter Kll) to represent the radiation illuminating 
the WA cloud. Assuming the plane parallel geometry, 
we ca rried the sim u lations usin g TI TAN p hotoionisation 
code (iDumont. Abrassart fc Collinl [2000 1. By assump¬ 
tion of constant total pressure, the matter of the WA 
is self consistently stratified without any additional re¬ 
quirements. The decrease of the radiation pressure is 
self consistently computed within photoionisation code 
TITAN and occurs due to absorption and emission i.e. all 
possible interactions of radiation field with matter. Since 
exact radiative transfer is computed, radiation pressure 
self consistently influences the matter structure. For the 
differences between escape probability mechanism and 
exact radiati ve transfer calculatio ns we refer readers to 
the paper by IDumont et al.l (120031 1. 

Since the radiation field influences the matter struc¬ 
ture and eventual thermal instabilities, S-curve, which 
is the dependence of the temperature T on the dynamic 
ionisation parameter S (Eqj2]for definition), can be fully 
constructed with the use of a single clo ud and with TI¬ 
TAN calculations (|Rbzanska et al.l 120061) . In this paper, 
we compare results obtained with two codes, titan and 
CLOUDY, for exactly the same physical input: spectral 
energy distribution (SED) of Mrk 509 as the incident 
radiation. The two codes differ in the way they treat 
the radiative transfer, and only with the use of titan 
we managed to compute thermally unstable regions and 
fully repr oduce deep minima in the observed AMD of 
Mrk 509 (iDetmers et al.l[2011f) . 

This result is a direct proof that thermal instabilities 
are responsible for the deep minimum in the column den¬ 
sities of AMD. Moreover, it provides strong evidence that 
the WA in AGN is compressed by radiation pressure, and 
it should be modelled with proper treatment of this ef¬ 
fect. 

The overall structure of the paper is as follows: sec¬ 
tion [2] contains the information about the incident ra¬ 
diation shape we used in photoionisation calculations. 
The description of the photoionisation codes and results 
of the photoionisation calculations are given in section [3] 
In section U we present the method of determining AMD 
from observations and calculating AMD from our mod¬ 
els. Finally, we compare our theoretical AMD with the 
observational data. The conclusions are presented in sec¬ 
tion [5] 

2. MRK 509 AND ITS SPECTRAL SHAPE 

Mrk 509 with redshift of 0.034397 dHuchra et al 119931) . 
is one of the best studied local AGN, with exceptionally 
high luminosity L(l-1000 Ryd)=3.2 x 10 45 erg s _1 clas¬ 
sified as a Syl galaxy. Usually it is considered to be one 
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Figure 1 . Spectral energy distribution of Mrk 509. The black 
line corresponds to the incident spectrum which is the input 
to the photoionisation codes titan and CLOUDY (see text for 
explanation). The red circles are the points adapted from Kll 
and normalised to incident flux at the corresponding inner 
radius of the cloud. 

of the closest QSO/Syl hybrids. iPeterson et al.l (|2004f ) 
have reported that Mrk 509 harbours a super massive 
black hole of mass 1.4 x 10 s M 0 . 

The WA in this source was extensively studied with 
600 ks RGS (reflection grating sp ectrometer) on board 
of XMM-Newton X-ray telescope (IDetmers et al.l 1201 It 
iKaastra et alJ 1201 2ft . The few tens of absorption lines 
from higly ionized metals were identified, allowing for 
determination of ionic column densities. Furthermore, 
equivalent hydrogen column densities were calculated for 
each ion, and absorption mea sure distribution was con¬ 
structed (IDetmers et al.|[20TTl) . Additionally, the upper 
limits on the location of different ionisation phases of the 
WA were estimated, using the variability method, which 
requires the assumption that any observ ed changes are 
cause d by changes in the ionisation states (IKaastra et al.l 
l20H . Those authors reported that observed ionisation 
components span the distance upper limits from 5 to 400 
pc. 

For any consistent photoionisation calculations of 
WA, the broad-band spectral energy distribution 

of the source is crucial. It_was shown by 

IRozahska. Kowalska fe Goncalvesl (120081 ) that if SED of 
an AGN is dominated by a soft component (i.e. an ac¬ 
cretion disk), the stability curve obtained during pho¬ 
toionisation calculations differs substantially from that 
obtained with the AGN SED dominated by a hard X- 
ray power-law component. The authors also explained 
that this is due to the fact that in the case of soft radi¬ 
ation entering the WA, bremsstrahlung is the dominant 
cooling in the cloud, while in the case of hard X-ray il¬ 
lumination, Compton scattering is responsible for high 
temperature equilibrium in the WA. A systematic study 
of how the stability curve depends on the shape of the 
incident SED was done bv iChakravortv et al.l (|2009f ). In 
particular, the authors discussed the influence of the soft 
X-ray excess on the shape of the stability curve in this 
paper. Additionally, chemi cal abundances affect the ion - 
isation balance of the WA dHcss. Kahn fe Paerelslll997l l , 
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and hence also the shape of the resulting stability curve. 
Nevertheless, for this paper we assume Solar chemical 
abundances, and we keep them constant during all pho¬ 
toionisation calculations. 

Recently, very detailed observed continuum spectral 
shape of Mrk 509 was published by Kll. The spectrum 
was obtained during their multi-wavelength campaign of 
several observations taken in various energy bands, with 
XMM Newton: RGS, EPIC (pn, MOS), and OM. Simul¬ 
taneously with XMM-Newton observations, data with 
INTEGRAL was also obtained to observe the hard X- 
rays. Additionally, they collected the data of Mrk 509 
with Chandra LETGS together with HST COS UV de¬ 
tection. Moreover, observation of Mrk 509 was done with 
Swift, using both the instruments: XRT and UVOT. For 
the final construction of SED, the authors used IR points 
from IRAS and Spitzer data, but it is unclear how much 
of IR emission is seen by the outflow. This is because 
most probably the wind is located closer to the black 
hole than the dusty torus from which most the IR radi¬ 
ation is originated. 

Kll allowed us to use those points for our project. The 
SED of Mrk 509 covers a wide range of wavelength band, 
which is essential for obtaining the ionisation balance 
needed for photoionisation modelling. All observed spec¬ 
tral points are presented in Fig.|T]as red circles. Mrk 509 
has one of the best spectra covering broad band among 
all the well studied AGN. Black line on the same figure 
demonstrates the incident SED used by us in photoion¬ 
isation calculations by titan and CLOUDY code. Both 
codes allow for the input of incident spectrum by points 
and the linear interpolation is done in order to fully de¬ 
fine incident continuum which goes into the simulations. 

3. PHOTOIONISATION CALCULATIONS 

Here, we study the detailed structure of the 
WA cloud using the pht oionis ation codes titan 
( Dumont. Abrassart fc CdollirJl200f)h and CLOUDY 13.02 
( Ferland et al.ll2013f) with the same initial physical condi¬ 
tions. Although both codes are used to study the clouds 
subjected to the incident radiation emanating from the 
central source, they do differ in their assumptions of ra¬ 
diative transfer treatments and number of spectroscopic 
lines considered. Nevertheless, the general idea is the 
same. Both codes solve radiation transfer through the 
gas in thermal and ionisation equilibrium with non-LTE 
equation of state. 

The most important difference between codes is that 
CLOUDY uses escape probability mechanism to solve ra¬ 
diative transfer, while TITAN uses Accelerated Lambd a 
Iteration (ALI) method (jGollin. Dumont fo Godctlf2004l l. 
which is more accurate in optically thick media. The 
detailed comparison of both radiative transfer meth- 
ods in the context o f X-ray absorbing gas was done in 
iDumont et alJ (|2003l) . Below, we show and compare the 
output of two codes in the case of our particular source 
Mrk 509 and its warm absorber. 

Following assumptions are made for this modelling in 
both photoionisation codes. The geometry of the cloud 
is set to the plane parallel. The expression for the ioni¬ 
sation parameter £, at the irradiated surface of the cloud 
relating to the distance from the illuminating source R 


is: 


where L lon is the luminosity of the ionising source and 
n is the hydrogen number density. The luminosity can 
be obtained from observations, whereas the number den¬ 
sity and ionisation parameter are free parameters in our 
calculations. 

3.1. TITAN computations 

tita n was developed by IDumont. Abrassart fc Gollinl 
( 2000 1 and was mainly designed to determine the struc¬ 
ture of temperature and ionisation state and the con¬ 
tinuum emission spectrum of a thick hot photoionised 
slab of gas. Using this code, the physical state of the 
gas is computed at each depth, assuming the local bal¬ 
ance between ionisation and recombination of ions, exci¬ 
tation and deexcitation, local energy balance and finally 
total energy balance, titan assumes the transfer of both 
continuum and lines using ALI method, which precisely 
computes line and continuum intensity self consistently. 
There are many options to control the physical proper¬ 
ties of the absorbing gas. Here, to calculate the stab ility 
curve which is compared to CLOUDY code (Sec. 13.311 . we 
use the assumption of constant density slab, but for the 
explanation of AMD in Mrk 509 we reject this assump¬ 
tion and compute more realistic constant total pressure 
cloud (see results in Sec.Q}. This exercise is made to en¬ 
sure the reader, that stability curves calculated by both 
codes agree for the same physical conditions of the gas. 

For the purpose of this paper, we make stronger as¬ 
sumption that the total pressure (i.e. gas and radiation 
pressure: P gas + P ra d) is held constant throughout the 
plane parallel slab of the gas, and the volume density 
and ionisation parameter across the cloud are stratified. 
The assumption of constant pressure allows the ionised 
gas to be naturally stratified due to illumination, as dis¬ 
cussed in (jRozanska et alJ l2006f> . The computations of 
all the structure are done iteratively until the conver¬ 
gence is reached. 

In the constant total pressure model, plane parallel ge¬ 
ometry is considered, and the radiation source is assumed 
to hit the cloud perpendicular to the surface. The param¬ 
eters of the model are the ionisation parameter £oj and 
the hydrogen number density no (both defined at the il¬ 
luminated cloud surface), the total column density Njj, 
and the incident SED. Although the ionisation parame¬ 
ter is well constrained in the surface of the illuminated 
cloud, ionisation properties are continuously changing as 
we go deeper into the cloud. This happens since all ther¬ 
modynamical parameters such as: gas temperature, gas 
density, and the radiation pressure do change with the 
depth of the cloud. To trace this behaviour, the dynamic 
ionisation parameter, 5 is defined as: 

—. S Lion 1 Pion Prad /.., 

47 xckT 47rcP 2 nkT cP &as P gas 

where Pi on is a flux affecting the cloud, c is a velocity of 
light and T is the temperature of the gas. Here, we do 
not follow the standard convention that P gas accounts 
only for hydrogen density number n. For fully ionised 
gas, with heavy element abundance, the total density 
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number equals ~ 2.3 n, and if the gas pressure is com¬ 
puted with total density, the factor 2.3 should be taken 
int o account, as it is in the basic definition of 5, given 
bv iKrolik. McKee fc Tartcrl (119811 ). Nevertheless, since 
the observers do not use the value of density while deter¬ 
mining AMD from observations, we compute the current 
value of £ using temperature, radiation pressure and the 
gas pressure, which directly comes out from our code. 

All thermodynamical parameters describing the cloud 
structure as temperature, radiation pressure, hydrogen 
and total density numbers are given as functions of the 
distance from illuminating source. Therefore, we can self 
consistently compute both ionisation parameter and its 
change with the distance from the cloud illuminated sur¬ 
face, using the value of radiation pressure computed from 
the second moment of radiation field. 

3.2. CLOUDY computations 

Cloudy is a well documented public code extensively 
used in the calculation of the photoionisation in many 
different clouds surrounding planetary nebulas, in the 
interstellar medium and Active Galaxies. It has an ex¬ 
tensive explanation of the physics of the plasma sub¬ 
jected to different perturbations given in the well organ¬ 
ised books: (Hazyl and Hazy2)0. The same as in titan 
runs for CLOUDY calculations, we assumed the open ge¬ 
ometry with thin plane parallel zones. All other input 
parameters as: the incident SED, the ionisation param¬ 
eter, the inner radius, the total column density of the 
slab and the hydrogen density at the surface of the illu¬ 
minated cloud are the same in both codes. We did this 
simulation for many grids of densities. The default op¬ 
tion in the CLOUDY is the constant density through the 
cloud. Nevertheless, there is an option to calculate cloud 
under constant pressure, and we check it and discuss in 
this paper. 

3.3. Stability curves for constant density clouds 

The general way for the analysis of the stability of the 
cloud subjected to the incident radiation from the cen¬ 
tral source is to study a stability curve, i.e. tempera¬ 
ture as a function of the ratio of the radiation pressure 
to the gas pressure, which is a dimensionless parameter 
S (Eq. El. Such defined stability curve can be computed 
in several ways. The most natural way is to compute 
a grid of constant density clouds located at the same 
distance from the illuminated source. By changing the 
cloud density, we change the ionisation parameter and 
so the gas temperature of the cloud. In Fig. [2j we com¬ 
pare stability curves computed in this way by two pho¬ 
toionisation codes CLOUDY (cyan triangles), and TITAN 
(magenta squares). Both codes use parameter selection 
of the clouds tabulated in Table [T] For each calculation, 
we considered the source luminosity 3.2 x 10 45 ergs s _1 
and the inner radius of the cloud, i.e., the distance be¬ 
tween the source and the illuminated face of the cloud is 
taken as log(r 0 /cm) = 17.25. 

In case of extended envelopes, the radiation pressure 
changes while we go deeper into the cloud due to geomet¬ 
rical dilution and absorption by ionised gas. Therefore, 
in principle, it is possible to calculate extended, spheri- 
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Figure 2. Equilibrium S-curves for constant density plane- 
parallel clouds obtained with the photoionisation codes titan 
magenta circles, and CLOUDY cyan triangles. Additionally, 
constant density single cloud with hh = 1 in spherical geom¬ 
etry computed with CLOUDY is presented by solid magenta 
line. 


Table 1 

Constant Density Cloud parameters. When the range of 
densities is shown that means the grid of clouds were 
calculated. The inner cloud radius in all cases is 
log(ro/cm) = 17.25. 


Codes Used 

710 (cm -3 ) 

£o (ergs cm s —x ) 

Geometry 

TITAN 

io 2 -io 12 

as L/nr z 

plane parallel 

CLOUDY 

10 2 -10 12 

as L/nr 2 

plane parallel 

CLOUDY 

1 

10 11 

Spherical 


cally symmetric cloud at constant density. While enter¬ 
ing deeper inside the cloud, the luminosity decreases as 
a distance square, and ionisation parameter of the con¬ 
secutive zone decreases. A full stability curve can be 
reconstructed also in this way, as presented in Fig. [2] by 
magenta solid line. 

All curves agree for temperatures down to 3 x 10 4 K. 
Below this temperature, atomic data which are different 
in the two codes may cause some inconsistencies. But 
since we concentrate here on the warm absorbers, we may 
consider that both codes work perfect in the considered 
temperature regime. 

3.4. Stability curves for constant pressure clouds 

As shown in Fig. [3] the constant pressure stability 
curves form titan and CLOUDY shows the similar trend 
except for the considerable difference in the thermal in¬ 
stability zone. Thermal instability zones are part of the 
stability curve where the slope is negative. The main 
difference between two codes here is that CLOUDY is not 
able to compute the structure at the zones of instability 
(there are no points on the unstable branch of stability 
curve). The stability curve computed by titan shows 
several points on unstable branch. Also there is a differ¬ 
ence in the range of ionisation parameter at which the 
instability occurs. Such difference may be due to the 
different approximations used in the code to solve ther¬ 
mal balance equation and the differences in the atomic 
database used. 
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Figure 3. Equilibrium S-curve obtained from CLOUDY (red) 
and titan (cyan). Both curves are produced with the number 
density no = 10 6 cm~ 3 and ionisation parameter £ = 10 6 erg 
cm s -1 at the irradiated surface of the cloud. The inner radius 
to the cloud is set to log(r D /cm) = 16.75 cm, to ensure that 
illuminating luminosity is equal to the observed one. The 
total column density (JVh) of the slab in both calculation is 
set to 1.5 x 10 23 cm -2 . 


the formula: 

N ion = A z@ j AMD/i on (log £)d(log £). (5) 

It is not possible to get the hydrogen column den¬ 
sity directly from observation since we do not see hy¬ 
drogen lines. It is derived using the ionic column den¬ 
sities of each ions assuming continuous distribution of 
the absorbers for se veral objects, including Mrk 509 
(jDetmers et all 12011 ). It should be noted that the hy¬ 
drogen column densities derived from the observations 
are nomin al a nd mo del dependent. For instance, in 
iDetmers et al.l (12011D paper, authors always used thin 
slabs, assuming that lines are unsaturated. Even adding 
the thin slabs together authors make a strong assump¬ 
tion that they are on the linear part of the curve of 
growth, and the fitted ionic column densities are always 
only lower limits. This kind of analysis is valid if the 
absorber is optically thin. For optically thick lines, the 
absorber can posses higher column density with the same 
strength of saturated lines. As a result the normalization 
of AMD fitted by thin slabs can be underestimated by 
one order of magnitude. Below, we present the modelled 
AMD, computed from the gas with continuous distribu¬ 
tion of ionisation parameters, as a result of the constant 
total pressure balance. 


4. ABSORPTION MEASURE DISTRIBUTION 

The distribution of the absorbing column in the line 
of sight is often described as Absorption Measure Dis¬ 
tribution (AMD). Almost always the large range or ion¬ 
isation state of gas located on the line of sight towards 
observer is seen in AGN with the warm absorber. A 
few tens of absorption lines is detected in high resolu¬ 
tion spectra of several AGN. Each line can be treated as 
an constant density absorber and using photoionisation 
code for the thin slab, we can calculate the ionisation 
parameter and ionic column density Won of the slab re- 
qu ired to produce su ch an obs erved absorbing feature. 
iHolczcr. Behar fc Kasnil 120071 ) formulated the expres¬ 
sion for quantifying the strength of absorption at these 
different ionisation states as: 


AMD 


dN H 

dflogO 


( 3 ) 


where, £ is defined by Eq. [TJ and Nh is the total hydro¬ 
gen column density along the line of sight. Furthermore, 
if solar abundances A Zq are assumed, the equivalent hy¬ 
drogen column density can be calculated separately for 
each ion using relation: 


N h 


N r 


./'ion (sin ax ) A Zq 


( 4 ) 


where /ion(£max) is the fractional ion abundance with 
respect to the total abundance of its element, and £ max 
is the maximum value of ionisation parameter required 
to produce the observed feature. In nature there is a 
distribution of gas with different £, which may account 
for absorption in the line, and in reality one must take 
into account the full dependence of /i on on £. This means 
that when ionic column densities for individual warm 
absorber are measured, the AMD can be derived from 


4.1. Modelled AMD for Mrk 509 

We used the titan photoionisation code to reproduce 
AMD for Mrk 509 assuming th at the absorbing clou d 
is in total pressure equilibrium (jRozanska et al.lI2006T) . 
When a cloud is irradiated with the incident SED of the 
source, the total pressure equilibrium requirement self 
consistently provides the distribution of the ionization 
parameters for which the AMD is calculated. Note, that 
when we compute AMD from simulations, we compute £ 
at each depth by relation: 

£ = 47i-cfcr-^- =4 ttc^-. (6) 

-Ugas(H) Tl 

Our model of constant pr essure reproduces th e AMD 
derived observationally by IDetmers et all ( 2011 ). Both 
discontinuities seen in the distribution of AMD in their 
paper are reproduced using our model. We made several 
exercises how the AMD behaves with different parame¬ 
ters and the results are as follows. 

For the detail investigation of the effect of the density 
on AMD structure, we made several runs with densities 
no = 10 5 , 10 6 , 10 8 and 10 9 cm -3 as presented in Fig. [4] 
(the density no = 10 7 cm~ 3 is removed from the figure 
for clarity). For lower values of the density at the illu¬ 
minated face of cloud no upto 10 8 cm -3 , the position of 
the dips remains the same . For higher densities the po¬ 
sition of dips in AMD cha nges. T his is con n ected wit h 
the fact reported by iRozariska. Kowalska fe Goncalvesi 
( 2008 ) that photoionisation models with high density 
similar to those in Broad Line Region (BLR) are not 
degenerated any more when SED of central source is 
dominated by strong optical/UV component. This is a 
case of Mrk509. The physical reason for this is that in 
the case of the hard incident spectrum all radiative pro¬ 
cesses (Comptonization, line heating/cooling) are linear 
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Figure 4. AMD structure for the same initial ionisation pa¬ 
rameter £o = 4.3 x 10 4 erg cm s~ 4 and different number den¬ 
sities at illuminated side of clouds marked in right bottom 
corner. 



Figure 5. AMD structure for the same initial number density 
no = 10 6 crrD 3 and different values of ionisation parameters 
at illuminated side of clouds marked in left upper corner. 

in density and the solution is determined by the ion¬ 
ization parameter, independently from the local density. 
However, for soft incident spectra the bremsstrahlung 
plays an important role, and the quadratic dependence 
of the bremsstrahlung (free-free) on the density breaks 
the degeneracy and introduces a dependence on the den¬ 
sity (|Rozariska. Kowalska fe Goncalvesll2008f) . 

As can be seen from the Fig. [2 when the densities 
are comparable to the Narrow Line Region (NLR) den¬ 
sities, the two dips remain prominent almost around the 
same range of ionisation states. As we move towards the 
regime where the density is comparable to the BLR den¬ 
sity, i.e. at no = 10 9 cm~ 3 , there is only one prominent 
dip present. Computations for higher densities show only 
one dip in AMD, but the amount of material in the ab¬ 
sorber produces much higher normalization of the mod¬ 
elled AMD than it is observed. It is difficult to say if the 
disappearence of this dip is connected with a switch in 
radiative cooling as it is valid in photoionisation models 


E 
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Figure 6. AMD structure for the same SED and different 
values of total column density considered. This is a case with 
no = 10 8 cm' 3 and £o = 3.4 x 10 3 erg cm s -1 . Absorbers with 
low column densities, less then 10 22 cm -2 , do not have enough 
matter to pass through two dips in AMD distribution. Those 
cases are indicated by blue diamonds and magenta squares. 


of different density. The line cooling may also be impor¬ 
tant. To confirm this result, the different SED should be 
considered and we plan to explore this problem in the 
future work. 

Overall normalization of the AMD does not change 
much for densities up to 10 8 cm -3 , but the normalization 
of dips is different. Dips are deeper for lower incident 
densities of the absorber. The depths of the two dips were 
the main criteria used in our comparison of modelled to 
the observed AMD presented below in Sec. 14.21 

In Fig. 0 we show the dependence of the AMD struc¬ 
ture on the initial ionisation parameter at the irradiated 
side of the absorber, £o- It can be seen that the location 
of the instability region remains at the same ionization 
range inside the absorber. The difference is only in the 
normalisation of the dips. 

There is no difference in absorption measure distribu¬ 
tion from the clouds of the same density affected by vari¬ 
ous luminosities if we assume the same ionization param¬ 
eter. This is obvious due to the fact that for the same 
product of the same energy flux L/R 2 is illuminat¬ 
ing the absorber. 

We also investigated the effect of different total col¬ 
umn densities of the cloud on AMD. It is shown in Fig. [G] 
and Fig. [7] for the number densities no = 10 s cm -3 and 
no = 10 9 cm -3 respectively at the illuminated side the of 
the absorber. As we can see from the Fig. 0 absorbers 
at lower total column density do not give two dips in 
AMD. Also their overall normalization is lower by about 
the half order of magnitude. If the total column den¬ 
sity is higher than 10 22 cm' 2 , the AMD has always two 
dips. The absorption measure distribution ends at lower 
ionization parameter for higher total column denisty of 
the constant pressure cloud. As seen in Fig. [7J when the 
column density is substantial, further increase of this pa¬ 
rameter does not change the distribution of AMD dips. 
This result is con sistent with the conclusion made by 
IStern et all (120141 ). 
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Figure 7. AMD structure for the same SED and different 
values of total column denisty considered. This is a case with 
no = 10 9 cm -3 and £o = 4.3 x 10 4 erg cm s -1 . 

4.2. Comparison with observations 

The direct comparison of an observed AMD with the 
modelled one is quite hard. This is because, for the given 
source we do observe only from several to several of tenth 
of absorption lines. Each line puts a point on AMD as 
described in the beginning of Sec.|4] Note, that in case of 
computed models, a few thousands of lines are included, 
therefore AMD has many more points. 

To date, observed AMD are found from individual 
observations of several objects. The general behaviour 
of AMD in case of five Syl galaxies was given by 
iBeharl (120091 1. All of those five galaxies: NGC 3783, 
NGC 5548, MCG-6-30-15, NGC 3516, NGC 7469, 
IRAS 13349+2438, reveal one dip in AMD located more 
or less in the same place between log£ = 0.8 — 2. The dip 
is not explained by RPC model presented in this paper. 
Additionally those objects are rather Syl galaxies with 
different SED than Mrk 509. The illuminating SED has 
the strong influence on AMD. For Sy 1.5 galaxy Mrk 509 , 
AMD reveals two dips as shown bv lDetmers et al.1 (1201 ill . 

The c omparison between the observed AMD for 
Mrk 509 (iDetmers et al.ll2011f) and our constant pressure 
model computed with realistic SED input taken from cur¬ 
rent multi-wavelenght studies (Kll) is shown in Fig'. [SI 
where observations are given by black histogram, and 
red triangles describe our model. Our calculations show 
the distribution of the ionization parameter much more 
densly. Furthermore, our column density presented in 
the figure by red triangles is the realistic column density, 
with radiation transmitted through all zones in the se¬ 
quence, instead of transmission through a co llection o f in¬ 
dependent zones, as in data modelling by (IDetmers et al.l 
1201 111 , as given above by Eq. [5] The lower panel of Fig. [5] 
shows absolute values of AMD normalization. It is obvi¬ 
ous that modelled normalization is by an order of magni¬ 
tude higher than the observed one in the stable zones. In 
the upper panel, we re-scale the normalization to show 
that dips are fully reconstructed by our model for the 
given incident density and ionization parameter. We ar¬ 
gue here, that the difference in the normalization occurs, 
since observed lines are saturated. Therefore, the derived 


ionic column densities from observations using thin slabs 
give only lower limits, and in reality can be one order of 
magnitude higher. To check our prediction we have to 
use saturated models to fit the data. They are not very 
common in the literature. 

There are few values of AMD normalization of several 
objects reported in the recent papers both from obser¬ 
vation and modelling. Summarizing, the agreement in 
the overall shape of the modelled and observed AMD 
dependence on the ionization parameter is very encour¬ 
aging but the normalization of the measured and mod¬ 
elled AMD remains as an issue. The observed AMD 
normalization of Mrk 509 corresponding to the ioniza- 
tion parameter £ = 8 x 10 2 erg cm s -1 is ~ 10 21 cm” 2 
(IDetmers et al.l [201 If) , less by a factor of ~ 30 than the 
AMD normalization from our best fitted constant total 
pressure model wi th TITAN (~ 3 x 10 22 cm -2 ). The 
AMD modeled by IStern et all (120141) within the frame 
of their RPC model were shown to be fairly universal, 
only weakly dependent on the model parameters, and the 
obtained normalization values were ~ 1.1 x 10 22 cm -2 . 
Measurements of AMD in six other objects (lBcharll2009f ) 
implied values of the order of ~ 4 x 10 21 cm -2 . Thus 
modelled values of AMD normalizations are frequently 
higher than the observed ones, and it is not clear whether 
the problem lies in the measurement method or in the 
modelling. 

Nevertheless, for the first time, the AMD dips can 
be shown in model computations. The discontinuities 
in AMD distribution are in the range of the ionisation 
parameter, log£ = 2 — 3 and log£ = 3 — 4. These 
ranges of £ corresponds to the ranges of temperatures 
T = (4 - 9) x 10 5 and T = (2 - 5) x 10 6 K, presented in 
Fig. [HI by the green solid line. This is a proof that the ab¬ 
sorbing material should be in pressure equilibrium and 
thin thermally unstable regions can be clearly visible. 
This comparison clearly shows that the warm absorber 
in Mrk 509 is the continuous cloud under constant total 
pressure. 

The two dipps are the result of the thermal instabil¬ 
ity operating in the irradiated medium. To show this 
effect we present the structure of main cloud parameters 
as the gas pressure, temperature, density and the ioni¬ 
sation parameter £ as the function of column density in 
Fig. [9j We show the plot of total cross section (a to t) 
versus column density in Fig. [TO] for the same case of 
cloud as in Fig. O Two cases of a sudden rise of the 
total cross section and a drop in the tempature can be 
clearly noticed, and we argue that they are responsible 
for the two dips in the observed AMD for Mrk 509. The 
discontinuity is not well resolved in the radiative trans¬ 
fer co des. As it was pointed out by iRozariska fe Czernvl 
(|1996f) thermally unstable zones are geometrically very 
thin and the correct temperature profile in the transition 
zone can be recovered only after including the electron 
conduction. Neither TITAN nor CLOUDY have this op¬ 
tion. However, the position of the discontinuity is well 
determin ed on the basis of rad iative transfer only, as was 
shown bv lCzernv et al.l (120091) and its extension depends 
on the model parameters. 

Our result puts strong constrains that the warm ab¬ 
sorber is build from the material of the density of the or¬ 
der of 10 8 cm -3 . For SED observed in Mrk 509, the value 
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Figure 8. The comparison of the simulated AMD structure - 
red triangles, with the observed one - black hi stogram. Data 
and e rror bars are taken from the paper by IDetmers et al.l 
(2012). The best model which roughly matches the data is 
computed for no — 10 8 cm -3 , £o = 3.4 x 10 3 erg cm s' 1 
and L = 3.2 x 10 45 erg s _1 . Bottom panel presents data with 
absolute normalization, while in the upper panel the modelled 
AMD is scaled down by a factor ~ 30 to match the observed 
points on the stable branch. The best fit model represents 
well the depth of the two dips. 

of density changes the nature of the AMD. For densities 
smaller than 10® cm -3 , we see two dips of different depth, 
as observed, and the best models compared to the data 
is for 10 8 cm' 3 . For densities of the order of 10 9 cm' 3 
one dip for higher ionisation parameter is substantially 
reduced. For the given SED we were unable to converge 
computations for hig her densities, b ut we plan to check it 
for objects from lStern et al.1 (l2014f) paper, with different 
shape of SEDs. Densities much lower than 10 8 lead to 
significant overprediction of the depth of f = 200 feature, 
as seen from Fig. [2 The discrepancy between the true 
modeled AMD normalization and the observed normal¬ 
ization, caused by the line saturation problems possibly 
broadens the acceptable density range. If indeed the sta¬ 
ble parts of AMD are underestimated in the data, but 
the dips are not, the true depth of the dips should be 
larger than in the data. Assuming the true depth of the 
low £ dip equals four decades, then the densities as low 
as 10 6 (but not lower) provide an acceptable solution for 
this part of AMD. However, the high £ dip is then not 



Figure 9. The structure of pressure (upper left), temperature 
(upper right), density number (bottom left), and ionisation 
parameter (bottom right), versus cloud column density for 
initial parameters no = 10 8 cm -3 , £o = 4.3 x 10 4 erg cm s' 1 
(solid lines) and no = 10 9 cm' 3 , £o = 4.3 x 10 4 erg cm s' 1 
(dashed lines). In both cases illuminating luminosity is L = 
3.2 x 10 4S erg s' 1 . 



Figure 10. Total cross section (<Ttot) relative to Thomson 
cross section (ot) as a function of cloud column density. The 
initial cloud parameters are as in Fig. [9] 

so well fitted. 

We note, that RPC model presented by iStern et all 
(|2014f ) should reproduce the observed AMD dips since 
thermal instability is clearly seen in constant pressure 
cloud presented in Fig. [3] by cyan triangles. It was shown 
that CLOUDY code fully repro duces therm al i nstab ilities 
under different circumsta nces (Rozahska & Czermll 1 9961 
iChakravortv et all 120121). In our opinion, the lack of 
AMD dip in the IStern et all (1201411 RPC model is due 
to the final binning of the AMD used by those authors. 
They have used final bins of 0.25 dex in £ which is 
too rough to show thin, unstable zones. Several very 
small di ps are ac tually present in the RPC model by 
IStern et al.l (|2014f ) (see their Fig. 2), but none of them 
explains observations. 

5. CONCLUSIONS 
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We performed photoionisation simulations of the warm 
absorber as a cloud illuminated by the energy emitted 
from the AGN centre. Computations were done assum¬ 
ing total pressure, P gas + P ra ,H to be constant, and with 
the use of TITAN code by Dumont. Abrassart fc Collin 

(eSsj. 

We computed absorption measure distribution for the 
case of Mrk 509, for which spectral energy distribu¬ 
tion is well known from the multi-wavelength long term 
monitoring. We found discontinuities in AMD in the 
range of the ionisation parameter,log^ = 2 — 3 and 
logf = 3 — 4, s i milar to that obtained observationally by 
iDetmers et al.l (120 Ilf) . Such observed discontinuity was 
often interpreted as the absence of o f ions i n th e ther - 
mally unstable regions dHolczer. Behar fc Kaspil 120071) , 
which can also be seen in the cooling curve in Fig. [3] 
These observed minimum is also described as two geo¬ 
metrically distinct regions along the line of sight repre¬ 
senting high ionisation r egion and low ionisation region 
jHolczei^^ehar^rKaspi 2007lh Contrary to the work by 
IDetmers et al.l ( 201 111 , instead of two discrete zones, we 
explain AMD as a continuous absorber under constant 
total pressure. 

Since the purpose of this paper is to reproduce the deep 
minima observed in the continuous distribution of the 
AMD, the modelled AMD is normalised to the observed 
one for convenience. The model requires higher values 
of column density to produce the structure all the way 
down to the low ionisation parameter. As a result of 
this, the normalisation of the modelled AMD is higher 
than the observed one. Nevertheless, the observed AMD 
are constructed with some approximations discussed in 
previous sections. 

The model presented in this paper clearly shows that 
only for the warm absorber with substantial density, we 
can reproduce AM D in Mrk 509. We point out here, that 
iStern et al.1 (I2014T) did not obtain AMD dips from their 
RPC model since their modeled AMD was computed 
with too strong binning. The number of dips observed in 
AMD may put constraints on the density of the absorb¬ 
ing material, which is a crucial parameter needed to in¬ 
dicate the wind location through the relation Eq. [T] Our 
work presents new density diagnostic by matching the 
model with exactly the same dips normalization to the 
observed AMD. The dips strongly reflect geometrically 
narrow thermally unstable regions. For higher densities, 
one dip disappears for a considered SED of Mrk 509. To 
confirm this result and to search for a parametrization 
of the shape of the AMD dips in terms of the incident 
n o, it is necessary to consider a sample of sources and a 
range of illuminating SEDs. We plan to investigate this 
issue in the future work. 
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